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Abstract 

In the context of the Higgs model involving gauge and Yukawa interactions with the 
spontaneous gauge symmetry breaking, we consider inflation with non-minimal grav¬ 
itational coupling, where the Higgs field is identified as inflaton. Since the inflaton quartic 
coupling is very small, once quantum corrections through the gauge and Yukawa interac¬ 
tions are taken into account, the inflaton effective potential most likely becomes unstable. 
In order to avoid this problem, we need to impose stability conditions on the effective 
inflaton potential, which lead to not only non-trivial relations amongst the particle mass 
spectrum of the model, but also correlations between the inflationary predictions and the 
mass spectrum. For concrete discussion, we investigate the minimal B — L extension of 
the Standard Model with identification of the B — L Higgs field as inflaton. The stability 
conditions for the inflaton effective potential fix the mass ratio amongst the B — L gauge 
boson, the right-handed neutrinos and the inflaton. This mass ratio also correlates with 
the inflationary predictions. In other words, if the B—L gauge boson and the right-handed 
neutrinos are discovered in future, their observed mass ratio provides constraints on the 
inflationary predictions. 
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1 Introduction 


Current understanding about the origin of our universe is that, for a very brief moment at 
the beginning, our universe went through a period of rapid accelerated expansion known as 
inflation. Inflation scenario [H IH |3l H] was originally proposed to solve serious problems in 
the Standard Big-Bang Cosmology, namely, the horizon, flatness and monopole problems. In 
addition and more importantly in the view point of the current cosmological observations, 
inflation provides a mechanism to create primordial density fluctuations of the early universe 
which seed the formation of large scale structure of the universe that we see today. In a 
simple inflation scenario, inflation is driven by a single scalar held (inhaton) that slowly rolls 
down to its potential minimum (slow-roll inhation). During the slow-roll era, the inhaton 
energy is dominated by a slowly varying potential, which causes the universe to undergo a 
phase of an accelerated expansion. Quantum huctuations of the inhaton held are stretched 
to macroscopic scales by inhation to yield the primordial density huctuations. After inhation, 
the inhaton decays to the Standard Model (SM) particles and the decay products heat up the 
universe (reheating). The success of big bang nucleosynthesis scenario requires the reheating 
temperature to be Tr > 1 MeV. 

Recently the Planck 2015 results [5] have set an upper bound on the tensor-to-scalar ratio 
as r < 0.11 while the best ht value for the spectral index (n^) is 0.9655 ± 0.0062 at 68% CL. 
Hence, the simple chaotic inhationary scenario with the inhaton potentials R oc and R oc 
are disfavored because their predictions for r are too large. Among many inhation models, Xcf)^ 
inhation with non-minimal gravitational coupling where 0 is inhaton, TZ is the scalar 

curvature, and is a dimensionless coupling) is a very simple model, which can satisfy the 
constraints by the Planck 2015 with ^ > 0.001 [6]. 

Given that we need interactions between SM particles and inhaton for a successful reheating 
of the universe, a more compelling inhation scenario would be where the inhaton held plays 
another important role in particle physics. As an example of such a scenario, we may consider 
the (general) Higgs model, where a scalar (Higgs) held plays the crucial role to spontaneously 
break the gauge symmetry of the model, and we identify the Higgs held as inhation. The SM 
Higgs inhation [3 Eli is nothing but this scenario, where the SM Higgs boson plays the role of 
inhaton with non-minimal gravitational coupling. Because of the observed Higgs boson mass 
of around 125 GeV, the SM Higgs ehective potential is likely unstable]^ If this is the cas^ 
the original SM Higgs inhation cannot work any more, and some extension is necessary [T2]rl 
However, we may apply the same idea to the general Higgs model and identify the Higgs held 
of the model (not the SM Higgs held) as inhaton in the presence of non-minimal gravitational 
coupling. For a simple example, see [T8] . 

As in the SM, the Higgs held in the general Higgs model has the gauge, Yukawa and quartic 
Higgs interactions. For complete analysis of inhation scenario in the Higgs model, we consider 

^ The stability of SM Higgs potential is very sensitive to the initial value of top quark pole mass, and we 
need more precise measurements for it [10]. See also m for the effect of possible Planck scale physics to the 
effective Higgs potential. 

Supersymmetric version of the Higgs inflation [H da uside] is free from this instability problem because 
of supersymmetry. It has been shown in that a large non-Gaussianity can be generated in this class of 
models. 
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the effective inflaton/Higgs potential by taking quantum corrections into account. In fact, we 
see that quantum corrections most likely cause an instability of the effective inflaton potential. 
Note that unless the non-minimal coupling (^) is very large, the quartic inflaton coupling is very 
small [6]. Hence, quantum corrections to the effective potential are dominated by the gauge and 
Yukawa interactions. We consider the renormalization group (RG) improved effective potential 
described as 

rw = Iaw <f\ (1.1) 

where 0 denotes inflaton, and A(0) is the running quartic coupling satisfying the (one-loop) 
RG equation of the form, 

~ Gg - Gy (1.2) 

cxln/i 

Here, g and Y are the gauge and Yukawa couplings, respectively, and Cg and Cy are positive 
coefficients whose actual values are calculable once the particle contents of the model are defined. 
Since the quartic coupling is very small, we have neglected terms proportional to A (A^ term 
and the anomalous dimension term). The solution to the RG equation is controlled by g and 
Y, which are much larger than A and independent of A. Therefore, we expect that unless 
the beta function is extremely small, the running inflaton quartic coupling A is driven to be 
negative in the vicinity of the inflation initial value, in other words, the effective potential has 
true minimum (far) away from the vacuum set by the Higgs potential at the tree-level. 

A simple way to avoid this problem is to require the beta function to vanish at the initial 
inflaton value (the stationary condition of A with respect to 0), namely, Cg g — Cy Y = 0. 
This condition leads to a relation between g and Y, equivalently, a mass relation between the 
gauge boson and fermion in the Higgs model. Since the Higgs quartic coupling at low energy is 
evaluated by solving the RG equation, in which the gauge and Yukawa couplings dominate, the 
resultant Higgs mass also has a relation to the gauge and fermion masses. The stability of the 
effective potential also requires the positivity of the second derivative of the potential, which 
leads to another constraint on the gauge and Yukawa couplings. In the slow-roll inflation, the 
inflationary predictions are determined by the slow-roll parameters defined with the potential 
and its derivatives, and therefore, the inflationary predictions have a correlation with the mass 
spectrum of the Higgs model. 

In order to explicitly show the mass relation and the correlation between the particle mass 
spectrum and inflationary predictions, we take the minimal B — L model as an example. This 
model is a very simple, well-motivated extension of the SM, where the global B — L (baryon 
number minus lepton number) in the SM is gauged. Three right-handed neutrinos and the B—L 
Higgs field (which is identified as inflaton) are introduced for the cancellation of the gauge and 
gravitational anomaly and the B — L gauge symmetry breaking, respectively. Associated with 
the B — L gauge symmetry breaking, the B — L gauge boson and the right-handed neutrinos 
acquire their masses. With the generation of the Majorana right-handed neutrino masses, the 
seesaw mechanism [19] for the light neutrino mass generation is automatically implemented in 
this model. Analyzing the RG evolutions of the B — L sector and the effective inflaton {B — L 
Higgs) potential, we show the particle mass spectrum and its correlation to the inflationary 
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predictions. Through the correlation, the Planck 2015 results provide us with constraints on 
the particle mass spectrum. 

This paper is organized as follows. In the next section, we briefly review the inflation 
with non-minimal gravitational coupling at the tree-level, and discuss the inflationary predic¬ 
tions in the light of the Planck 2015 results. In Sec. 3, we introduce the minimal B — L extension 
of the SM and calculate the RG improved effective Higgs potential. We show the particle mass 
spectrum derived from the stability conditions of the effective potential and its correlation to 
the inflationary predictions. We then compare our results for various values of the non-minimal 
gravitational coupling ^ with the Planck 2015 results. In Sec. 5, we discuss reheating scenario 
in the B — L Higgs inflation for the completion of our inflationary scenario. Sec. 6 is devoted 
to conclusions. 


2 Non-minimal Xcf)^ inflation at tree-level 


In the Jordan frame, the action of our inflation model is given by (hereafter we always work in 
the Planck unit, Mp = Mpi/\f^ = 1, where Mpi = 1.22 x 10^® GeV is the Planck mass) 


Sj = I 




(2.1) 


where /{cj)) 
potential is 


(1 -1- ^0^) with ^ being a positive, dimensionless parameter, and the inflaton 


V{<t>) = 



( 2 . 2 ) 


Using the conformal transformation, gEfiu = f{(p)gtiu, the action in the Einstein frame is descried 
as 


Se = 


(L^Xy/-gE 



T^e + 2 eO'Y 


PW)]. ■ 


(2.3) 


In the Einstein frame with a canonical gravity sector, we describe the theory with a new inflaton 
held (cr) which has a canonical kinetic term. The relation between a and the original inhaton 
held (j) is given by 




l + (l + 60e</>^ 
(i+e02)2 


(2.4) 


where a prime denotes the derivative with respect to 0. 

The inhationary slow-roll parameters in terms of the original scalar held (0) are expressed 
as 


e(0) 

viP 

C(0) 


U vyV 

2 [VEa'J ’ 

VeP'Y VEp'f ’ 

( vk \ ( yp o vk^'" \ 

\VEa') VVE(a')3 VE{a'Y VE{a'Y VE{a'Y) ’ 


(2.5) 
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Figure 1: The inflationary predictions for iV = 50 (dot-dashed) and iV = 60 (solid) in the non- 
minimal inflation. The top panels show r vs. ^ (left) and vs. ^ (right). The bottom-left 
panel shows A vs. The inflationary predictions Ug and r for various values of ^ are depicted 
in the bottom-right panel, along with the contours at the confidence levels of 68% and 95% 
given by the results of Planck 2015 (solid) and Planck-|-BICEP2/Keck Array (dot-dashed) |5]. 


where Ve is the potential in the Einstein frame in terms of original held 0 given by 


Ve(c^(0)) 


1 A 0^ 

4(l + e02)2 • 


The amplitude of the curvature perturbation A-ji is given by 


( 2 . 6 ) 



Ve 

2A7r'^e 


ko 


(2.7) 


which should satisfy = 2.195 x 10 ® from the Planck 2015 results [5] with the pivot scale 
chosen at ko = 0.002 Mpc“^. The number of e-folds is given by 


2 h. ye(0) \d4>/ ’ 


( 2 . 8 ) 


where 0/ is the inflaton value at horizon exit corresponding to the scale ko, and 0e is the 
inflaton value at the end of inflation, which is defined by max[e(0e), |r7(0e)|] = 1. The value 
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of N depends logarithmically on the energy scale during inflation as well as on the reheating 
temperature, and is typically around 50-60. 

The slow-roll approximation is valid as long as the conditions e <C 1, jr/l 1 and C 1 
hold. In this case, the inflationary predictions, the scalar spectral index Ug, the tensor-to-scalar 
ratio r, and the running of the spectral index a = are given by 

rig = I — Qe + 2r], r = 16e, a = 16er] — 24e^ — 2^. (2.9) 

Here the slow-roll parameters are evaluated at ^ = 0/. 

Fig. [T] shows inflationary predictions for the non-minimal inflation at the tree-level for 
= 50 (dashed-dotted) and iV = 60 (solid). Top panels show r vs. ^ (left) and Ug vs. ^ 
(right). Both r and Ug show asymptotic behavior for both small and large ^ values. In the 
minimal inflation limit with ^ = 0, we obtain r ~ 0.31 (0.26) and Ug ~ 0.942 (0.951) for 
= 50 (60). The plots also show that for a larger e-holding number, we obtain a larger Ug 
while a smaller r. The bottom-left panel shows the tree level quartic coupling A as a function of 
Note that A is very small unless ^ 1. The inflationary predictions for Ug and r for various 
values of ^ are depicted in the bottom-right panel along with the results from the measurements 
by Planck 2015 and Planck-|-BICEP2/Keck Array [5]. We see that the inflationary predictions 
for ^ > 0.001 are consistent with the observations. 

3 Running B — L Higgs inflation and stability of inflaton 
potential 



SU(3), 

SU(2)i 

U(1)k 

U(1)b l 


3 

2 

+1/6 

+1/3 

Ur 

3 

1 

+2/3 

+1/3 


3 

1 

-1/3 

+1/3 

4 

1 

2 

-1/2 

-1 

AT?* 

1 

1 

0 

-1 


1 

1 

-1 

-1 

H 

1 

2 

-1/2 

0 

L 

1 

1 

0 

+2 


Table 1: Particle contents of the minimal B—L model. In addition to the SM particle contents, 
the right-handed neutrino AT?* {i = 1,2,3 denotes the generation index) and a complex scalar 
ip are introduced. 

In order to investigate the Higgs inflation with the stabilized inflaton potential, in this paper 
we take the minimal B — L extension of the SM as an example, where the anomaly-free U(1 )b_l 
gauge symmetry is introduced along with a scaler field ip and three right-handed neutrinos AT?*. 
The particle contents of our model are listed in Table 1. This model requires three generations 
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of right-handed neutrinos to cancel all the gauge and gravitational anomalies. The B — L gauge 
symmetry is broken by the vacuum expectation value (VEV) of (f in its Higgs potential of 

H(|<p|) = A (<^V-• (3.1) 

Associated with the gauge symmetry breaking, the right-handed neutrinos acquire their Majo- 
rana masses through the Yukawa interaction, 

1 3 _ 

£ D ^ Y^NR^^NR + h.c., (3.2) 

i=\ 


where we have taken the degenerate mass spectrum for the right-handed neutrinos, for sim¬ 
plicity. After the B — L symmetry breaking with the Higgs VEV (ip) = vbl/V^, the particle 
masses are given by 

mz' = 2 g vbl, riY^R = ^ H vbl, vbl- (3.3) 

v2 


Let us now consider the B — L Higgs inflation scenario. The action in the Jordan frame is 
given by 


^tree 



-t/(l4=l)K + - V(M) 


(3.4) 


where D^ = d^ — i2gZ'^^ and /(|<^|) = 1 -l- 2^(pt(p, and (p = {vbl + 0)/\/2 in the unitary gauge 
with the physical Higgs hied 0 identihed as inhaton. In the Einstein frame, the RG improved 
effective inhaton potential at the one-loop level is given by izni 

VeW = iA(4) (3.5) 


where <I> = + ^0^, and we have neglected vbl which is much smaller than the Planck 

mass. The RG equations of the couplings at the one-loop level are given by [18] 


= i(4s + 32)y. 

16>rV^ = -6g^Y + Y^ 

a/i 2 

167rV^ = + 2) A^ - (48^2 - 6Y^)X + 96^ - 3Y^ (3.6) 


where a factor s dehned as 


(1 -|- 


l + (l + 60eh" 

is assigned to each term in the RG equations associated with only the physical Higgs boson 
loop corrections [T] |8]@ We hnd that this s-factor has no ehect in our numerical analysis for a 
parameter region satisfying g^, Y^ 3> A, and one may hx s = 1 as a good approximation. 


® There are a few different prescriptions for computing quantum corrections in the presence of the non- 
minimal gravitational coupling [^. For recent, detailed computations of quantum corrections, see [20] and their 
results of 1-loop beta functions with the s-factor. 
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Figure 2: The RG evolution of the inflaton quartic coupling A for various values of g and R, 
with the hxed value of .^ = 1. The initial value of A(<F/) is hxed by the tree-level analysis. In the 
left panel, the solid, the dashed and the dot-dashed lines show the running quartic couplings for 
g = 0.01, 0.011 and 0.009, respectively, with the hxed value oiY = 0.0237. In the right panel, 
the solid, the dashed and the dot-dashed lines show the running couplings for Y = 0.0237, 
0.0214 and 0.0261, respectively, with the hxed value oi g = 0.01. The stability condition of 
= 0 is satished for g = 0.01 and Y = 0.0237, and the corresponding running couplings 
are depicted as the solid lines. 

Let us now investigate the stability of the ehective inhaton potential. In our analysis 
throughout this paper, we set the initial values of A to be the one obtained in the tree-level 
analysis at the initial inhaton value (p = (pi, equivalently, = (pi/■yjl + ^cpj. Then, we consider 
the RG improved ehective inhaton potential by taking into account the RG evolution of the 
quartic coupling with the initial condition at (pi. As we have seen in the previous section, the 
inhaton quartic coupling is very small unless 3> 1. Hence the beta function of the quartic 
coupling is approximately given by 

/9a = [{I8s^ + 2) A^ - (48/ - 6Y^)X + 96/ - 3^^] ~ ^ (96/ - 3^") , (3.8) 

when g'^, Y"^ :$> A0 The RG evolution is controlled by g and Y, which are independent of the 
initial value of the inhaton quartic coupling. Fig. [2] shows the RG evolution of the inhaton 
quartic coupling in the vicinity of the initial inhaton value for various values of g and Y with 
a hxed .^ = 1. In the left panel, the solid, the dashed and the dot-dashed lines denote the 
RG evolutions for g = 0.01, 0.011 and 0.009, respectively, with the hxed value oi Y = 0.0237. 
In the right panel, the solid, the dashed and the dot-dashed lines denote the RG evolutions 
for Y = 0.0237, 0.0214 and 0.0261, respectively, with the hxed value oi g = 0.01. Note that 
/5 a( 4>/) = 0 is satished with the parameter choice for the solid lines, g = 0.01 and Y = 0.0237. 
We can see from Fig. |2]that if the condition of = 0 is violated even with ±10 % deviations 
for the values oi g ot Y, the running quartic coupling quickly becomes negative in the vicinity 

® In this paper, we are interested in this case, otherwise the beta function is so small that the inflaton quartic 
coupling is almost RG invariant. Although the tree-level analysis is valid in this case, the gauge and Yukawa 
couplings are too small to yield any impacts in the experimental point of view. 
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of $/ (see the dashed and the dot-dashed lines). This fact indicates that the B — L gauge 
symmetry breaking vacuum at 0 = vbl is unstable and the effective potential develops a true 
vacuum with a negative cosmological constant. The quantum corrections through the gauge 
and Yukawa coupling completely change our inflationary scenario from the one at the tree-level. 




Figure 3: The RG improved effective potential for various values of g and a hxed .^ = 0.1. In 
the left panel, the effective potentials for g = 0.041 (solid) and 0.046 (dot-dashed) are shown. 
For g = 0.046 > gmax = 0.0425, the local minimal is developed. The vertical dashed line 
indicates <F = 4)/ (in the Planck unit). In the right panel, the effective potentials for various 
values of gmin < 9 < gmax are shown. 

In order to avoid this instability, we impose not only the condition of = 0 but also 
dPx/d^ > 0 at <F/. From the hrst condition, the Yukawa coupling Y is determined by the 
gauge coupling, which we take as a free parameter in our analysis, along with the others, ^ and 
Vbl- The second condition ensures that the effective potential is monotonically increasing in 
the vicinity of $/, and yields a lower bound on > gmin- When we analyze the global structure 
of the effective potential, we can notice that there exists an upper bound on g < gmax- For a 
large g > gmax, the effective potential develops a local minimum at <F < $/, so that the inflaton 
held will be trapped in this minimum after inhation. A second inhation then takes place until 
the vacuum transition from this local minimum to the true B — L symmetry breaking vacuum. 
To avoid this problem, the parameter region is restricted to be in the range of gmin < 9 < Qmax- 
Fig. [3]shows the effective potential for various values of the gauge coupling (g) for hxed ^ = 0.1. 
In the left panel, the solid line depicts the ehective potential for g = 0.041, while the dot-dashed 
line is for g = 0.046. In this example, we hnd the upper bound as gmax = 0.0425. We can see 
that the ehective potential develops a local minimum for g = 0.046 > gmax- For various values 
of g < gmax, the ehective potentials are shown in the right panel. 

4 Inflationary predictions and low energy observables 

Under the stability conditions, gmin < 9 < gmax with various value of we now calculate 
the inhationary predictions with the ehective inhaton potential. Since we refer the results in 
the tree-level analysis for A($/) for a hxed ^ and impose the stability condition /3 a($/) = 0, 





our prediction for the tensor-to-scalar ratio r is the same as the one obtained in the tree-level 
analysis. However, the RG evolution of the inflaton quartic coupling alters the other inflationary 
predictions, Ug and a, from those obtained in the tree-level analysis, because they are calculated 
by the second and third derivatives of the effective potential (see Eqs. fl2.5p and fl2.9p ). 

Let us hrst derive a analytic formula for the deviations of Ug from the tree-level prediction. 
In the effective inflaton potential, the inflation quartic coupling is not a constant, but a function 
of <h = With the stability condition = 0, we calculate the derivatives of 

the effective potential as 




A'($) 


<t>=4>i 4 


+ (i*)' . . = Ll*) (i*)' 


(j)=4>i 4 


+iA'(<h) (<h^)' + ^A($) (<h^)' 

4>=4>i 2 4‘=4>i 4 


4>=4>i 


= V'e 


tree 


aa"($) 

4 




+ V^ 


tree 


(4.1) 

(4.2) 


tree 


and Ve 


are evaluated to be 


tree 


where a prime denotes the derivative with respect to 0, 

the same as those in the tree-level analysis, and we have used X'{(f>i) oc /SA(‘h/) = 0 under the 
stability condition. Rewriting d/dcf) in terms of d/d^, we have 


6 ( 0 /) = etree, = Vtree + 


(4.3) 


where etree and rjtree are the slow-roll parameters evaluated in the tree-level analysis, and 


(l-e4>^)^ A(4>,) /)a(4>7) 

l + A($7) l + 4>/A($/)‘ 


(4.4) 


Here, a dot is the derivative with respect to <h. Since is non-zero, rj is deviated from its 
tree-lvel value rjtree- Using Eqs. fl3.6p and (13.81) . we obtain 

167r^$//3A(4>/) ^ 4608 / + 72 - 30 U® ~ 768 (^9 - 5^2^ (4.5) 

where in the last expression, we have used the stability condition ~ at $/. Finally, we 
arrive at an expression for the spectral index as (see Eq. (12.91) )!^ 


ng = nT" + 2^rj ~ nf"" -f 


6(9 -5\/2) 


TT^ 


{1 + A($7)’ 


(4.6) 


For a hxed <h/ and A(<h/) are determined by the tree-level analysis, so that the inflationary 
prediction is controlled by the gauge coupling g, which is the free parameter in our analysis. 

In our model, there are only three free parameters, vbl and g, with a hxed e-folding 
number N = 50/60. Once we hx ^ and N, <h/ and A(<h/) are hxed by the tree-level analysis for 

^ In the same way, we can express the running of the spectral index, a, in terms of $/, A and g. However, 
as we will see in the following numeral analysis, the predicted a values are found to be very small and always 
consistent with the Planck 2015 results. Thus, we omit the expression for a. 
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inflation, and the inflationary predictions except for r are controlled by the gange coupling g 
with its relation to the Yukawa coupling Y led by the stability condition = 0. In the B — 

L model, the particle mass spectrum is determined by the gauge, Yukawa and inflation quartic 
couplings at the scale vbl (see Eq. fl3.3l) h which are obtained by solving the RG equations in 
Eq. fl3.6p from /i = $/ to p = vbl- Since A($/) is very small, its RG evolution is determined by 
g and Y as we can see from its RG equation. Gonsidering all of these facts, we expect that there 
exist a nontrivial mass relation in the particle spectrum and a nontrivial correlation between 
the inflationary predictions and the particle mass spectrum. In the following, the results of our 
numerical analysis will show such non-trivial relations. For simplicity, we hx vbl so as to yield 
mz' = ‘2g{vBL)vBL = 3 TeV, to be consistent with the current results from the search for Z' 
boson resonance at the Large Hadron Collider [21], [22] . 



Figure 4: Inflationary predictions for various values of ^ with the inputs of g varied in the 
range of gmin < 9 < Qmax, along with the contours given by the experiments (same as in Fig. [1]). 
The horizontal solid lines from top to bottom correspond to the results for ^ = 1.5 x 10“^, 
2.1 X 10“^, 2.9 X 10“^, 4.1 X 10“^, 6 xl0“^, 0.01, 0.02, 0.1 and 1. The diagonal dashed line 
denotes the inflationary predictions in the tree-level analysis. 

Fig. m shows the resultant inflationary predictions for a variety of ^ values with the input 
values of g in the range of gmin < 9 < 9max at {N = 60), along with the contours given by 
the Planck 2015 results. As we have discussed above, the prediction for the tensor-to-scalar 
ratio is the same as the one in the tree-level analysis, while the predicted spectral index is 
altered by quantum corrections. In Fig. 0] we can see that for < 1, the results show sizable 
deviations for g ~ gmax from those at the tree-level analysis depicted as the diagonal dashed 
line. Interestingly, the Planck 2015 results provide upper bounds on g, which are more severe 
than gmax for ^ < 0.001. In our numerical analysis, we can see that d*/ approaches l/\/^ from 
a smaller value as we increase ^ > 1, and hence the deviation of the predicted Ug value from 
the one in the tree-level analysis becomes smaller as we can see from Eq. fl4.6p with the limit 
1 ^ 0 . 

In order to see the inflationary predictions as a function of g, we show our results in Fig.[5]for 
^ = 0.0029, 0.01 and 0.1, with N = 60. The top-left and top-right panels show the inflationary 
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Figure 5: Inflationary predictions and low energy observables for various values of 5 ^ < gmax 
and ^ ~ 0.0029, 0.01 and 0.1 (corresponding r ~ 0.108, 0.045 and 0.008), respectively. We have 
fixed iV = 60 and mz' to be 3 TeV. Qmax increases with increasing Top two panels show the 
inflationary predictions for spectral index Ug (left) and the running of spectral index dns/d\nk 
(right), for decreasing r (top to bottom). Bottom left panel shows mass of the inflaton (top 
to bottom) for decreasing r. Bottom right panel shows mass ratio mNn/mz' (top to bottom) 
for increasing r values. 


predictions of Ug and a = dng/dlnk as a function of g in the range of gmin < g < gmax- 
The prediction for r is the same in the tree-level analysis, and r ~ 0.108,0.045 and 0.008, 
respectively, for ^ = 0.0029, 0.01 and 0.1. For a larger value of g-max becomes larger. As g 
is lowered, the predicted Ug value approaches the tree-level prediction. By numerically solving 
the RG equations for the couplings in Eq. fl3.6p for a fixed g, we obtain the particle mass 
spectrum with niz' = 3 TeV. The mass ratio rriN/mz' is shown in the bottom-left panel, while 
the bottom-right panel shows the inflaton mass. In both panels, the solid lines from left to right 
correspond to the results for ^ = 0.0029, 0.01 and 0.1, respectively. The resultant mass ratio 
is almost independent of but shows a splitting for g ^ 5 x 10“^. For such a very small g, its 
corresponding Y determined by the stability condition is also very small, and hence both g and 
Y are almost RG invariant and the mass ratio is determined by Y{^j)/g{^i). However, in this 
case, the condition g^yY"^ 3> A is no longer valid, and Y determined by /9A(*h/) = 0 depends 
on the input A values. This is the reason why the bottom-left panel shows the splitting among 
three solid lines for g ^ 5x 10“'^. For g > 0.005, the RG evolution of A is mainly determined by 
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g and Y in its beta function, since A($/) is extremely small. For a very small g value < 0.005, 
the effect of g and Y on the RG evolution of A becomes negligible, and X{vbl) — A($/). Since 
we have hxed mz' = ‘^g{vBL)vBL = 3 TeV, vbl = 1-5 TeY/gigjBh) — 1-5 TeV/ 5 f(<h/), and the 
inflaton mass becomes larger proportionally to l/g{^i) as shown in the bottom-right panel. 
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Figure 6: Same as Fig. [5] with ^ = 0.1, but for two different values of = 50 (dashed lines) 
and = 60 (solid lines). 

Same as Fig. [5] but for A^ = 50 and 60 with ^ = 0.1 is depicted in Fig. [6l The dashed lines 
denote the results for N = 50, while the solid lines for N = 60. The inflationary predictions 
show a sizable deference for the two different N values, as shown in the tree-level analysis in 
Sec. 2. On the other hand, the particle mass spectrum weakly depends on N . 

We show the results for the mass spectrum for large ^ values and A^ = 60 in Fig. [7] The 
solid lines from left to right corresponds to the results for ^ = 1, 10, 50, 150, 500 and 1000. 
For 9 ^ 0.05, the resultant mass ratio in the left panel shows ^-dependence. This is because 
S> A is no longer valid for such a small g value, and the Y value determined by 
depends on A(<h/). We also show the inflation mass spectrum in the right panel, which show a 
similar behavior as the results in the bottom-right panel in Fig. [5l As we have seen in Fig. 01 
the inflationary predictions for > 1 are close to those obtained in the tree-level analysis, 
rig — 0.968 and r ~ 0.003. 

Finally, we show in Fig. [H]the contour plots (diagonal solid lines) for the inflationary pre¬ 
dictions with various fixed g{vBL) values, along with the results shown in Fig. 01 Here, for a 
hxed g{vBL) value, we calculate the inhationary predictions from the ehective inhaton poten¬ 
tial for various values of We have used mz' = 3 TeV The diagonal solid lines correspond to 
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Figure 7: Mass spectrum for various value of ^ = 1,10, 50,150, 500 and 1000 from left to right. 
The left panel shows the mass ratio rriNji/mz', while the right panel shows the mass of the 
inflaton. Here we have hxed mz' = 3 TeV. 

gi^BL) = 0.0184, 0.0216, 0.023, 0.026, and 0.0425 from left to right. 



n, 


Figure 8: Inflationary predictions for various hxed g{vBL) values along with the results shown 
in Fig. m Here, for a hxed givsi) value, the inhationary predictions are calculated for various 
values of ^ taken in Fig. HI We have used mz' = 3 TeV. The diagonal solid lines correspond to 
g[vBL) = 0.0184, 0.0216, 0.023, 0.026, and 0.0425. 


5 Reheating after inflation 

Any successful inhation scenario requires the transition to the Standard Big Bang Cosmology 
after inhation. This happens via the decay of inhaton into the SM particles during the era 
inhaton is oscillating around it potential minimum, and the decay products then reheat the 

universe. We estimate the reheating temperature after inhation Tr by using V = H = 
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where T is the inflaton decay width, and is the effective degrees of freedom for relativistic 
SM particles when the reheating occurs, so that 


T„ ~ 0,2 


(5.1) 


From the success of big bang nucleosynthesis, we take a model-independent lower bound on the 
reheating temperature as Tr > 1 MeV. 

As we discussed in the previous section, the inflaton is much lighter than the Z' boson and 
right-handed neutrinos. Thus, its decay width to the SM particles through off-shell process 
mediated by the heavy particles is too small to satisfy the lower bound Tr > 1 MeV. An 
efficient reheating process is possible when the inflaton in general has a coupling with the SM 
Higgs doublet {H) such as 

C D -A' [H^H) ((pV) • (5.2) 

Although this coping is crucial for the reheating process, we assume A' <C 1 not to change our 
analysis for the B — L model in the previous sections. Since the inflaton quartic coupling is 
very small and hence the inflaton is light, it is most likely that the inflaton can decay to the 
SM particles only through its mixing with the SM Higgs boson through the A' coupling. After 
the B — L and electroweak symmetry breakings, we diagonalize the scalar mass mass matrix of 
the form, 

'h 

where h is the SM Higgs boson, and 01 and 02 are the mass eigenstates. The relations among 
the mass parameters and the mixing angle are the following: 


cos 6 sin 9 
— sin 6 cos 6 


(5.3) 


2vblVsm>^’ = {ml - ml) tan26>. 


=^h- {ml - ml) 
ml^ =ml+ {ml - ml) 


sin^ 9 

l-2sin^e’ 
sin^ 8 

l-2sin^e ■ 


(5.4) 


where vsm = 246 GeV is the Higgs doublet VEV, and mh and m^ are the masses for h and 0, 
respectively. 

The left panel in Fig. |9] shows the mass eigenvalue m^^ as a function of sin^ for various 
m^p values with mh = 125 GeV. In the 6 = 0 limit, mip^ = m^, while m^ = mh for 6 = 7r/2. 
Although in this plot we show the results by using Eq. (15.dh . we only consider the case with 
6* 1 as mentioned above, otherwise our results obtained in the previous sections are changed 

in the presence of a sizable A'. For 6* 1, the mass eigenstate 02 (0i) is almost identical to 0 

(h). 

The inflaton can decay to the SM particles through the mixing with the SM Higgs boson. 
We calculate the inflaton decay width as 




(6.5) 


where Thimp^) is the SM Higgs boson decay width if the SM Higgs boson mass were The 
reheating temperature is then evaluated by Eq. fl5.ip . For various inputs of m^, we show the 
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Figure 9: Left: the mass eigenvalue as a function of sin6' for various values with 
nih = 125 GeV. We have used Eq. fl5.4p . Right: the reheating temperature as a function of 
sin 6* for various values of = 0.32, 3.0 and 130 GeV (solid lines from bottom to top), along 
with the dotted horizontal line for the lower bound of Tr = 1 MeV. 


resultant reheating temperature in the right panel of Fig. [HI The solid lines from bottom to top 
denote the results for = 0.32, 3.0 and 130 GeV, respectively, along with the dotted horizontal 
line for the lower bound of Tr = 1 MeV. The sharp drop for each solid line corresponds to the 
fact that m <^2 becomes zero for a certain sin0 valne as shown in the left panel. We find that 
the universe is snfficiently heated up for a mixing angle in the snitable range shown in Fig. jH] 
Finally, we check the theoretical consistency of our analysis. When we introdnce the conpling 
in Eq. fl5.2p . the beta fnnction of the inflaton qnartic conpling is modified to 

levr^/?;, ^ Ibvr^/?;, + 2A'2. (5.6) 


In order not to change our results in the previous sections by the introduction of A', A'^ should 
be negligibly small in the beta fnnction. We then impose a condition 


A'2 


< 1 . 


(5.7) 


For m^, we obtain from Eq. fl5.4p 

(5,8) 

VblVsM mz'VSM 

where we have used mz' = 3 TeV, rrih = 125 GeV and vsm = 246 GeV. Gombining this 
eqnation with Eq. fl5.7p . we obtain 


^>6xlO“^0. (5.9) 

From Figs. |5l [7| and HI we can see that this condition is satisfied for a large potion of the 
parameters space. 
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6 Conclusions 


Inflationary universe is the standard paradigm in modern cosmology, which not only solves 
the problems in the Standard Big Bang Cosmology, but also provide the primordial density 
fluctuations necessary for generating the large scale structure of the present universe. As a 
simple and successful inflationary scenario, we have considered the inflation with non- 
minimal gravitational coupling. With a suitable strength of the non-minimal coupling, the 
inflationary predictions of this scenario becomes perfectly consistent with the Planck 2015 
results. 

It is more interesting if the inflaton can also play some crucial role in particle physics. 
We have considered the general Higgs model with the gauge and Yukawa interactions with 
the spontaneous gauge symmetry breaking. In the presence of the non-minimal gravitational 
coupling, the Higgs field can also play the role of inflaton. The analysis with the Higgs po¬ 
tential at the tree-level leads to the inflationary predictions consistent with the cosmological 
observations. However, once we take quantum corrections, the effective inflaton potential most 
likely becomes unstable. This is because the inflaton quartic coupling is extremely small in 
a large portion of the parameters space and the effective potential is controlled by the gauge 
and Yukawa couplings independently of the quartic coupling. In the renormalization group 
improved effective potential, we see that the running quartic coupling becomes negative in the 
vicinity of the initial inflaton value, indicating the instability of the effective potential. In order 
to avoid this problem, we have imposed the stability condition of vanishing the beta function of 
the inflation quartic coupling. This condition leads to a non-trivial relation between the gauge 
and fermion masses. Since the renormalization group evolution of the inflaton quartic coupling 
is mainly controlled by the gauge and Yukawa coupling, the inflation mass at low energy is 
determined by the couplings. Therefore, the mass spectrum of the gauge boson, fermion and 
inflation shows a non-trivial relation. 

Since the inflaton potential is modihed from the tree-level one, the inflationary predictions 
are altered from those obtained by the tree-level analysis. Although the prediction of the 
tensor-to-scalar ratio remains the same under the condition of the vanishing beta function, the 
predictions for the scalar spectral index and the running of the spectral index can be significantly 
altered. The fact that the effective potential is controlled by the gauge and Yukawa couplings 
implies a correlation between the inflationary predictions and the particle mass spectrum. 
Therefore, the observables at the gauge symmetry breaking scale correlate with the inflationary 
predictions which determined by physics at an extremely high energy compared to the gauge 
symmetry breaking scale. 

By taking the minimal B — L extension of the Standard Model as a simple example, we 
have shown such a non-trivial relation in the particle mass spectrum driven by the stability 
condition of the effective inflaton potential. We also have calculated the inflationary predictions 
from the effective potential and found their dependence of the B — L gauge coupling. Therefore, 
the new particle mass spectrum of the B — L model, once observed, has an implication to the 
inflationary predictions. On the other hand, more precise measurements of the inflationary 
predictions yield a constraint on the B — L particle mass spectrum. 

For completeness, we have also investigated reheating after inflation. Since the inflation 
is lighter than the Z' boson and the right-handed neutrinos, its reheating process through 
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the heavy particles are not efficient, and the resultant reheating temperature is too low to 
be consistent with the bound from big bang nucleosynthesis. We then introduce a coupling 
between the inflaton and the Standard Model Higgs doublet. Through the mixing with the 
Standard Model Higgs boson, the inflaton can decay into the Standard Model particles and the 
universe can be heated up with a sufficiently high reheating temperature. We have found that 
this happens with a sufficiently small coupling between the inflaton and the Higgs doublet and 
such a small has essentially no effect on our analysis for the particle mass spectrum and the 
inflationary predictions. 
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